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Report  No.  HE-150-108 


HEAT  TRANSFER  FROM  SPHERES  TO  A RAREFIED  GAS  IN  SUBSONIC  FLOB 


1.0  INTRODUCTION 

Convective  heat  transfer  from  spheres  to  a rarefied  gas  has  been 
treated  analytically  and  experimentally  in  Refs.  1,  2,  3 and  4.  From  the 
experimental  results  of  Ref.  4,  the  predicted  effect  of  gas  rarefaction 
on  the  heat  transfer  in  the  slip  flow  region  was  in  part  substantiated,  in 
that  the  heat  transfer  coefficient  decreased  below  its  continuum  value, 
with  increasing  rarefaction.  This  effect  is  apparently  (Ref.  1)  primarily 
due  to  the  temperature  jump  condition  at  the  surface,  analogous  to  the  slip 
velocity  at  the  surface,  which  occurs  at  low  pressures.  This  temperature 
jtwip  introduces  an  effective  thermal  contact  resistance  at  the  surface. 


The  experimental  results  of  Ref.  4 were  determined  under  super- 
sonic flow  conditions.  Since  there  existed  a curved  detached  shock  wave, 
some  doubt  remained  as  to  the  effect  of  the  Bhock  wave  and  possibly  other 
compressibility  phenomena  on  the  results.  It  was  considered  desirable  to 
obtain  heat  transfer  data  from  spheres  at  subsonic  Each  numbers  in  the  slip 
flow  region. 


It  is  the  purpose  of  this  report  to  present  overall  average  heat 
transfer  data  from  spheres  to  a rarefied  air  stream  for  a range  of  vari- 
ables, 0.1414  0.69  and  1.74Re4l24>  and  to  present  a simplified  analysis 
which  predicts  the  trends  of  this  data  by  a rarefaction  correction  to  the 
continuum  solution.  An  expression  for  the  non-dimensional  overall  average 
heat  transfer  coefficient  for  spheres  in  subsonic  flow  is  obtained  in  the 
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The  flow  regions  investigated  here  and  in  Ref.  4 are  shown  on  Fig.  1 in 
comparison  with  the  range  of  earlier  investigations  reported  in  Ref.  5. 

2.0  ANALYSIS 

If  one  considers  the  temperature  Jtaip  in  the  slip  flow  region 
as  an  effective  thermal  contact  resistance  at  the  surface  over  and  above 


- 2 - 


the  thermal  resistance  due  to  the  viscous  boundary  layer,  then  the  heat 
transfer  coefficient  at  low  pressures  can  be  determined  to  the  first  order 
by  a correction  to  the  ccntisuvsn  heat  transfer  coefficient  at  the  same 
Reynolds  number. 

Consider  a surface  placed  in  a continuum  flow  tie  that  he&t  is 
transfered  between  the  surface  and  the  stream.  At  any  particular  instant 
the  surface  temperature  is  assumed  uniform  and  equal  to  t 8 . The  heat 
flux  will  fcs 


a_± 

d y 


y = o 


(2.1) 


where  is  the  gradient  of  temperature  norma?  to  the  surface,  K is  the 
heat  conductivity,  hc  is  the  convective  heat  transfer  coefficient,  and  te 
is  the  equilibrium  temperature  of  the  surface.  The  subscript  y-o  iden- 
tifies the  properties  of  the  air  evaluated  at  the  surface  and  zero  super- 
script signifies  continuum  conditions. 


Now  if  the  air  stream  is  rarefied  while  holding  the  Reynolds  num- 
ber constant,  the  additional  resistance  to  heat  transfer  due  to  the  tempera- 
ture Jump  will  present  itself  so  that  we  may  define  a total  convective  heat 
transfer  coefficient  incorporating  both  the  thermal  boundary  layer  and  tem- 
perature ;ump  resistances  by  writing 

^ = nc  '^"w  ~ ^e)  (2o2) 

where  and  t are  the  surface  temperature  and  equilibria  surface  tem- 
perature respectively,. 


The  first  order  temperature  Junp  boundary  condition  may  be 

written  as 
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where  t 
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and 
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are  the  tempereture 
ture  in  the  layer  of  gas  immediately  adjacent  to 
tity  X is  the  temperature  juap  distance  defined 
(Ref.  6)  by  the  formula 
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and  gradient  of  tempera- 
the  surface.  The  quan- 
from  kinetic  theory 


(2.3) 


or: 


(2.4) 
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with  X M le48  , cL  the  aeccoitodation  coefficient,  f the  ratio  of 

■ s * 

specific  heats,  Pr  the  Prandtl  number,  and  0 = 1.996  -~h-  . 

CX  f + \ 

. c 

This  first  order  temperature  jump  condition  requires  thatTe*fe  . 

It  is  now  further  asst»ied  that  q , the  heat  transfsrcd  between  the  sphere 

at  temperature  and  the  gas  at  temperature  1\  for  the  case  of  a surface 

temperature  jump,  is  the  same  as  the  heat  transfered  between  a sphere  at 

temperature  _ and  the  gas  at  temperature  T,  under  conditions  of  no 

temperature  jxsap.  In  other  words*  Eq.  2.1  is  assumed  valid  vith  t =t  , 

W j - o 

Combinations  of  Eqs.  2.1,  2,2  and  2.3  give 


(2.5) 


Substituting  Eq.  2.4  in  2.5  yields 
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(2.7) 


■where  Mu  and  Mu  are  -the  local  Husselts  numbers  for  contlnmm  and  rarefied 
flows  respectively. 

Eq„  2,7  indicates  that  experimental  data  could  be  correlated  by  a 

*lot  of  VtJG  - nTjO  Ter4US  -dV'  ’ “e  ?ig- A- 

3.0  gPERTMKBTAL  KqTTTFireNT  AND  PROCEDURE 

3.1  The  required  low  density  gas  stream  was  provided  by  the  No.  3 

Wind  Tunnel  located  at  Berkeley,  The  constructional  features  and  gen- 
eral operating  characteristics  of  the  tunnel  are  contained  in  F.ef.  7. 
An  axi-symmetric,  subsonic  nozzle  was  employed  in  the  investigation, 
giving  a range  of  Mach  ntmbera,  0,1  to  0.69,  and  a static  pressure 
range  of  36  to  3300  microns  Hg.  The  nozzle  used  had  an  exit  diaaeter 
of  9 inches.  The  design  of  this  nozzle  iB  described  in  Ref.  8. 

The  four  spheres  and  the  impact  pressure  probe  were  counted  on  a 
rotary  selector  (see  Appendix  C,  Ref,  8)  which  in  turn  was  supported 
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in  the  tunnel  teat  chamber  by  a traversing  mechanism  (Ref.  7)  capable 
of  axial,  lateral,  and  vertical  movement  with  respect  to  the  nozzle. 
With  this  versatile  mounting  system,  any  of  the  five  objects  on  the 
rotaiy  selector  could  be  moved  into  the  flow  field  as  required,  with- 
out alteration  of  the  flow. 

3.2  Experimental  Method 

The  overall  average  heat  transfer  coefficients  were  determined  by 
the  transient  technique  described  in  Ref.  4=  Briefly,  the  method  is 
as  follows:  The  rate  of  change  of  heat  in  a small  body  of  sensibly 

uniform  internal  temperature  distribution  is  equated  to  the  heat  less 
by  convection  and  radiation  to  give 

- C wV  = hc  A (j  - fe  ; + nr  A (t  - ta) 


with  the  boundary  conditions 

i = t*  'tr  = o 

t -+e'  r « 


where  t,  is  the  initial  sphere  temperature  and  t ' is  the  equilibrium 

I C 

temperature  for  the  combined  effects  of  radiation  and  convection.  If 
the  difference  between  the  absolute  surface  temperature  and  the  abso- 
lute temperature  of  the  surroundings  is  not  great,  Eq.  3.1  may  be  in- 
tegrated for  h„  and  h not  functions  of  the  temperature  to  yield, 

W I 
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where  the  equilibrium  temperature  f. 


t 

is, 
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nr 
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and  the  equilibrium  temperature  due  to  convection  heat  transfer  alone 
is  given  from  Eq.  3.3  as 


It  is  necessary  to  determine  the  radiation  heat  transfer  coeffi- 
cient hf  independently  by  making  transient  cooling  runs  at  no  flow 
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with  tunnel  pressures  reduced  tc  the  order  of  0.1  am  Eg,  so  that  free 
convection  currents  are  eliminated.  The  heat  balance  in  this  case 
becomes 
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Making  use  of  the  temperature  time  histories  on  a semi-logar- 
ithmic plot  as  suggested  by  Eqe.  3.2  and  3.6,  the  radiation  heat 
transfer  coefficient,  hr  , and  the  sum  of  the  radiation  and  con- 
vection coefficients  (h  + h ) are  obtained.  The  quantity  t is 
measured  ty  means  of  a thermocouple  placed  on  the  tunnel  wall  and 
f 1 is  that  equilibrium  temperature  attained  by  the  sphere  after 
a very  long  time.  The  equilibrium  temperature  for  convection  only 
may  then  be  calculated  from  Eq.  3.4.  The  heat  transfer  coefficient 
due  to  convection  alone,  hc  , is  the  arithmetic  difference  between 
the  two  experimentally  determined  coefficients  (jn,Q+^r')  and  h • 

3.3  Spheres 

The  spheres  were  made  of  silver  and  were  mounted  on  hollcar  drawn 
glass  stings  as  shown  in  Fig.  2.  B and  S No.  40  iron  and  constantan 
wires  were  made  into  thermocouples  and  soft  soldered  into  the  centers 
of  the  spheres  and  were  then  led  out  through  the  glass  stings  to 
larger  sires  on  the  sting  support  leading  to  the  recording  potenti- 
ometer. The  spheres  are  further  described  in  Ref.  4. 

3.4  Experimental  Procedure 

The  sphere  in  question  was  heated  initially  to  a temperature  of 
approximately  120°-150°F  by  means  of  a radiation  furnace  consisting 
of  a IX  watt  lamp  encased  in  several  radiation  shields  and  mounted 
near  the  nozzle  in  the  test  chamber.  This  radiant  energy  was  suffi- 
cient to  give  the  sphere  a 10*  temperature  rise  in  approximately  10 
seconde.  The  sphere  was  then  rsmoved  from  the  furnace  and  traversed 
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into  the  air  stream  on  the  nozzle  centerline  with  its  forward  stagnation 
point  tangent  to  the  exit  plane  of  the  nozzle.  The  temperature-time 
history  of  the  sphere  for  this  position  was  recorded  by  means  of  the 
recording  potentiometer. 

3,5  Instrumentation 

All  pressures  were  measured  by  a precision  U-tube  manometer  (Kef,  9) 
to  an  accuracy  of  - 1 microns  Hg  for  pressures  up  to  400  microns  and 
1 1/4  per  cent  for  pressures  above  400  microns. 

Sphere  temperatures  were  measured  by  an  iron-constantan  thermo- 
couple and  recorded  by  a Leeds -Norturup  potentiometer  to  an  accuracy 
of  - .017  millivolts  - 0.5*?  and  2.5  psr  cent  in  time.  The  tunnel  wall 
temperature,  tQ  , and  settling  chamber  temperature,  tQ  , were  measured 
by  eopper-constantan  thermocouples  located  in  the  settling  chamber  and 
recorded  with  a Brown  16  point  recording  potentiometer  to  - .03  milli- 
volts (-  1 .4*? ) . 


In  addition  to  the  above  special  instrumentation,  the  standard 
tunnel  equipment  and  instrumentation  as  described  in  Ref.  7 was  utilized, 

4.0  REDUCTION  OF  EXPERIMENTAL  DATA 


4.1  Flow  System 


The  flow  parameters  required  in  the  data  presentation  are  liach 
number  and  Reynolds  numbers  evaluated  at  the  center  of  the  nozzle  exit 
plane.  These  values  were  determined  in  the  following  manner:  The 

static  pressure,  , was  measured  by  a wall  orifice  in  the  constant 
area  section  of  the  nozzle  and  the  total  or  impact  pressure  pj  was 
measured  by  a Bource-shaped  impact  probe  (0.300"  O.D.)  at  the  nozzle 
exit  plane.  Agreement  between  this  impact  pressure  and  the  reservoir 
pressure  allowed  the  use  of  the  following  isentropic  formula  for  the 
determination  of  Mach  number: 


jLll 


for  air  with  Y'  * 1,4 


r- 1 


(4.1) 
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a.  236 
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(4.2) 
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The  stagnation  temperature  was  assumed  identical  to  the  measured 
settling  chamber  temperature.  As  a result,  the  state  temperature  was 
determined  from  the  Mach  number  by  the  adiabatic  formula 


and  for  Y*  * 1,4 


Assuming  that  perfect  gas  relationship  holds, the  above  formulas 
for  M and  » together  with  the  measured  value  of  and  the  coeffi- 
cient of  viscosity  for  dry  air  at  temperature  (Ref.  10)  enabled  the 
Reynolds  number  to  be  calculated  by  the  formula 


Rex  --  6.64  X 10 
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In  Table  I are  tabulated  the  tunnel  flow  conditions  for  the  various 
rims.  No  no37le  blocking  corrections  were  applied. 


4*2  Thermal  System 


The  reduction  of  the  thermal  data  was  accomplished  by  means  of 
the  method  outlined  in  Section  3.2  to  result  in  the  values  of  the 

convection  coefficients  h_  . The  final  results  are  shown  in  Table  II 

C m hcD 

in  the  dimensionless  form  of  the  Nusselts  number,  !\]U  = r-4*—  . wherein 

' av  ke 

the  thermal  conductivity  of  the  air  is  evaluated  at  the  equilibrium 
temperature  of  the  sphere.  A typical  plot  of  log  (f  - versus  t“ 
appears  in  Fig,  3. 


5.0  SOURCES  CF  ERROR 


5.1  Flow  System 

Errors  incurred  in  the  determination  of  the  Mach  and  Reynolds 
numbers  were  primarily  due  to  the  inaccuracy  of  the  instrumentation 
used  to  measure  the  pressures  and  temperature  required  in  Eqs,  4.2 
and  4,5=  The  masiBum  relative  errors  of  those  quantities  vary  from 
0.1  to  25  per  cent  and  are  listed  with  the  quantities  themselves  in 
Table  I. 

5.2  Thermal  System 


(4.3) 

(4.4) 

(4.5) 


The  errors  attendant  to  the  determination  of  the  heat  transfer 
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poefficients  reside  in  the  measurement  of  the  t im e -temp era ture 
history  of  the  cooling  body  and  in  the  separate  determination  of  the 
equilibrium  temperature  tg  1 „ These  errors  are  thus  determined  by 
the  least  count  of  the  recording  potentiometers , The  largest  least 
count,  - .52®?,  when  associated  with  the  smallest  temperature  dif- 
ference measured  and  used  in  the  computation,  5®R,  gives  a possible 
error  of  10  per  cent,,  In  addition,  the  error  due  to  time  measurement 
is  2.5  per  cent.  The  radiation-conduction  correction  which  is  ac- 
counted for  experimentally  in  independent  runs  is  subject  to  the  tsamw 
errors  due  tc  instrumentation.  Since  the  radiation-conduction  cor- 
rection may  amount  to  50  per  cent  of  the  total  heat  transfer  coef- 
ficiont(bc+  ) in  the  worst  case, the  maximum  uncertainty  in  the 
determination  of  h^,  may  possibly  be  25  per  cent.  This  uncertainty 
in  hc  will  be  transmitted  directly  to  the  Nusselts  number  with  any 
other  error  incurred  in  the  value  of  the  thersaal  conductivity  of  the 
air. 


6.0  DISCUSSION  OF  EXPERIMENTAL  RESULTS 


The  analysis  given  in  Section  2.0  results  in  ar,  expression  for 
the  dimensionless  local  heat  transfer  coefficient  for  the  case  of  slip 
flow  by  virtue  of  an  effective  thermal  contact  resistance  due  tc  the  tem- 
perature jump  boundary  condition.  The  result  appears  as  a function  of  the 
dimensionless  local  heat  transfer  coefficient  for  continuum  flow  and  the 
local  Hach,  Reynolds  and  Prandtl  numbers.  It  is  given  as 

Nj  u - 


or 
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The  experiments!  data  is  in  the  form  of  overall  average  heat  transfer  co- 
efficients on  spheres  and  exact  prediction  of  these  results  cannot  be  ex- 
pected from  the  foregoing  simplified  analysis.  However,  one  could  expect 
that  So.  6„1&  would  give  a method  of  correlating  the  data,  insofar  as  it 

V A 

reduces  to  the  continuum  value  for  small  p0  V and  gives  a solution  for 

large  nr-  vr  which  resembles  the  free  molecule  solution  of  Ref.  3,  i.0=, 
k i ■ j RePrp/v^ 

5o~  - > for  S“-ailM> 
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For  our  particular  application  to  overall,  average  heat  transfer 
from  spheres,  we  shall  utilize  the  solution  from  Ref.  2 as  the  specifica- 
tion for  which  gives 


N u " = 
cv 


_QY 


k i 


n = 


(6.2) 


where  =\:,pr  Re  Pr 
VJ,,v  the  average  velocity  over  the  sphere . The  experimental  data  is  cor- 
related  as  ^ — j^Tr  versus  and  appears  in  Fig.  4 with  ri  = 1.3. 
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The  value  of  9 = 1.3  was  chosen  merely  to  match  the  solution  to  the  data 


for  small 


IS/! 


The  data  correlates  as  a straight  line  of  slope  equal 


Re  Pr 

to  3»42  as  determined  by  the  method  of  least  squares.  It  may  be  seen  in 
Fig.  4 that  batter  than  90  per  cent  of  all  the  experimental  data  lies  within 
a 10  per  cent  error  by  Mu  of  this  line.  This  is,  in  general,  within  the 
accuracy  of  measurement  of  the  data. 


The  above  empirically  determined  constant,  3.42,  enables  or.e  to 
write  an  expression  for  the  heat  transfer  for  spheres  in  a slightly  rare- 
fied subsonic  air  stream  as 


N U 


a* 


I + 3.42 


The  superimposed  plot  of  Eq.  6.3  as  well  as  the  theoretical  results  for 
spheres  from  Refs.  2 and  3 are  shown  for  comparison  with  the  experimental 
data  in  Fig.  5 as  Nusselts  number  versus  Reynolds  number  with  liach  number 
sb  the  parameter.  This  figure  presents  graphically  the  growing  influence 


of  the  rarefaction  parameter 


M 


Re 


upon  the  heat  transfer  and  also  illus- 


trates the  ability  of  the  proposed  rarefaction  correction  to  match  the 
continuum  and  free  molecule  flow  solutions. 


(6.3) 


In  Fig.  6 analytical  results  of  the  continuum  low  Reynolds  num- 
ber and  free  molecule  flow  regions  are  shown  for  comparison  with  known 
existing  experimental  data  for  heat  transfer  from  spheres. 

Temperature  recovery  factors  have  not  been  presented  because 
the  instrumentation  used  prevented  obtaining  an  acceptable  accuracy  in 
the  measurement  of  this  quantity. 
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7.0  RESULTS 

7.1  Experimental  average  overall  heat  transfer  coefficients  were 
obtained  for  spheres  in  a rarefied  subsonic  air  stream  where  the  Mach 
mmber  and  Reynolds  number  varied  from  0.1  to  0.69  and  1.7  to  124, 
respectively. 

7.2  A semi-empirical  formulation  of  the  Nusselts  number  for  spheres 
in  a rarefied  subsonic  air  stream  is  obtained  by  correcting  the  con- 
tinuum solution  for  an  effective  thermal  contact  resistance  due  to 
the  temperature  jump  boundary  condition. 

7.3  The  presented  data  follow  the  same  trends  toward  the  free  mole- 
cule solution  as  was  first  shown  in  Ref.  4 for  spheres  in  rarefied 
supersonic  flows. 
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TABLE  I 


SUBSONIC  NOZZLE  FLOW  CONDITIONS 


Tunnel 
Run  No. 

i 

Flow  Rate 

p. 

(miorons 

Hg) 

^0 

(°R) 

if 

m 

Re 

(per  inch) 

(per 

Error 

ceno) 

( lbs /hr ) 

u 

Re 

236  (9 ) 

20 

193 

538 

.69 

110 

+-  .1 

- .6  ! 

236(b) 

20 

374 

539 

.37 

116 

- 1.3 

- 1.6 

237 

20 

363 

535 

.37 

106 

- 1.3 

- 1.6 

238 

2.3 

36 

543 

.59 

17 

- 5.7 

- 8.5 

239 

2.3 

136 

537 

.17 

17 

*17 

-17 

240 

2.3 

63 

543 

.35 

17 

- 7,7 

- 9.3 

241 

2.3 

62 

538 

.37 

18 

- 7.7 

- 9.3 

242(a) 

2.3 

61 

536 

.37 

18 

- 7.7 

- 9.3 

242(b) 

2.3 

241 

536 

.10 

18 

-25 

-25 

243 

2.3 

236 

538 

.10 

17 

-25 

-25 

244 

50 

3293 

539 

.10 

248 

-16 

*-16 

245 

40 

1305 

543 

.21 

205 

- 3.6 

- 3.9 

246(a) 

40 

738 

538 

.37 

207 

- 1.2 

- 1.5 

24.6(b) 

40 

440 

543 

.59 

211 

- .5 

t .8 

247 

■■  - 

30 

570 

530 

.37 

166 

- 1.2 

- 1.5 

TABLE  II 
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FIG. 3-TYPICAL  PLOT  OF  SPHERE 
TEMPERATURE  WITH  TIME 
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FIG.  5-AVERAGE  HEAT  TRANSFER  COEFFICIENTS  FOR  SPHERES 
IN  SUBSONIC  FLOW  TO  RAREFIED  AIR 


